Controlled Synchronization of One Class of Nonlinear Systems under Information 
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Output feedback controlled synchronization problems for a class of nonlinear unstable systems 
under information constraints imposed by limited capacity of the communication channel are an- 
alyzed. A binary time-varying coder-decoder scheme is described and a theoretical analysis for 
multi-dimensional master-slave systems represented in Lurie form (linear part plus nonlinearity de- 
pending only on measurable outputs) is provided. An output feedback control law is proposed based 
on the Passification Theorem. It is shown that the synchronization error exponentially tends to zero 
for sufficiantly high transmission rate (channel capacity). The results obtained for synchronization 
problem can be extended to tracking problems in a straightforward manner, if the reference signal 
is described by an external (exogenious) state space model. 

The results are applied to controlled synchronization of two chaotic Chua systems via a commu- 
nication channel with limited capacity. 
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I. INTRODUCTION 

Analysis and control of the behavior of complex in- 
terconnected systems and networks has attracted con- 
siderable recent interest. The available results signif- 
icantly depend on models of interconnection between 
nodes. In some works the interconnections are modeled 
as delay elements. However, the spatial separation be- 
tween nodes means that modeling connections via com- 
munication channels with limited capacity is more real- 
istic. 

Recently the limitations of control under constraints 
imposed by a finite capacity information channel have 
been investigated in detail in the control literature, see 
[E HHH> SHi Q and the references therein. It has been 
shown that stabilization of linear systems under informa- 
tion constraints is possible if and only if the capacity of 
the information channel exceeds the entropy production 
of the system at the equilibrium (Data Rate Theorem) 
[1 0, HI]. In [1, [j| a general statement was proposed, 
claiming that the difference between the entropies of the 
open loop and the closed loop systems cannot exceed the 
information introduced by the controller, including the 
transmission rate of the information channel. 

For nonlinear systems only a few results are available 
in the literature |l El El El EE El- In the above 
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papers only the problems of stabilization to a point are 
considered. 

In the control literature there is a strong interest in 
control of oscillations, particularly in controlled synchro- 
nization problems EH, El El EH- However, results of 
the previous works on control systems analysis under in- 
formation constraints do not apply to synchronization 
systems since in a synchronization problem trajectories 
in the phase space converge to a set (a manifold) rather 
than to a point, i.e. the problem cannot be reduced to 
simple stabilization. Moreover, the Data Rate Theorem 
is difficult to extend to nonlinear systems. 

The first results on synchronization under information 
constraints were presented in [T{|, where the so called 
observer-based synchronization scheme [H [2lJ was con- 
sidered. In this paper we extend the results of [19( and 
analyze an output feedback controlled synchronization 
scheme for two nonlinear systems. A major difficulty 
with the controlled synchronization problem arises be- 
cause the coupling is implemented in a restricted man- 
ner via the control signal which is computed based on 
a measurable innovation (error) signal which has been 
transmitted over a communication channel. Key tools 
used to solve the problem are quadratic Lyapunov func- 
tions and passification methods [H l23j |. To minimize 
technicalities we restrict our analysis to Lurie systems 
(linear part plus nonlinearity depending only on measur- 
able outputs). 

The paper is organized as follows. The controlled syn- 
chronization problem is described in Section |TTJ The cod- 
ing procedure used in the paper, is presented in Sec. IIHI 
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The main results are presented in Section IIVI where an 
exponential convergence of the synchronization error to 
zero is established. An example showing synchroniza- 
tion of the chaotic Chua systems is presented in Section 
IVl Final remarks are given in the Conclusion. Auxiliary 
results are presented in the Appendices. 



II. DESCRIPTION OF CONTROLLED 
SYNCHRONIZATION SCHEME 

Consider two identical dynamical systems modeled in 
Lurie form (i.e. the right hand sides are split into a linear 
part and a nonlinear part which depends only on the mea- 
surable outputs). Let one of the systems be controlled by 
a scalar control function u(t) whose action is restricted by 
a vector of control efficiencies B. The controlled system 
model is as follows: 

x{t) = Ax(t) + B<p{ yi ), yi (t) = Cx(t), (1) 
z(t) = Az{t) + B^(y 2 ) + Bu, y 2 (t) = Cz(t), (2) 

where x(t), z(t) are rt-dimensional (column) vectors of 
state variables; yi(t), y 2 {t) are scalar output variables; 
A is an (n x n)-matrix; B is n x 1 (column) matrix; C 
is an 1 x n (row) matrix, ip(y) is a continuous nonlin- 
earity, acting in the span of control; vectors x, z stand 
for time-derivatives of x(t), z(t) respectively. System 
(TIj) is called the master (leader) system, while the con- 
trolled system ^ is called the slave (follower) system. 
Our goal is to evaluate limitations imposed on the syn- 
chronization precision by limiting the transmission rate 
between the systems. The intermediate problem is to 
find a control function IA(-) depending on the measurable 
variables such that the synchronization error e(t), where 
e(t) = x(t) — z(t) becomes small as t becomes large. We 
are also interested in the value of the output synchro- 
nization error eit) — yi(t) — 2/2 (0 = Celt). 

A key difficulty arises because the error signal between 
the master system and the slave systems is not available 
directly but only through a communication channel with 
a limited capacity. This means that the synchronization 
error e(t) must be coded at the transmitter side and code- 
words then transmitted with only a finite number of sym- 
bols per second thus introducing error. We assume that 
the observed signal e(t)) is coded with symbols from a fi- 
nite alphabet at discrete sampling time instants tk = kT s , 
k = 0, 1, 2, ... , where T s is the sampling time. Let the 
coded symbol e[k] = S(tk) be transmitted over a digital 
communication channel with a finite capacity. To sim- 
plify the analysis, we assume that the observations are 
not corrupted by observation noise; transmission delay 
and transmission channel distortions may be neglected. 
Therefore, the discrete communication channel with sam- 
pling period T s is considered, but it is assumed that the 
coded symbols are available at the receiver side at the 
same sampling instant tk = kT s , as they are generated by 
the coder. Assume that zero-order extrapolation is used 



to convert the digital sequence s[k] to the continuous- 
time input of the controller e(t), namely, that e(t) — e[k] 
as kT s < t < (k + 1)T S . Then the transmission error is 
defined as follows: 

S £ (t) = e(t) - e(t). (3) 

On the receiver side the signal is decoded introducing 
additional error and the controller can use only the signal 
e(t) = e(t) — S e (t) instead of e(t). 

We restrict consideration to simple control functions 
in the form of static linear feedback 

u(t) = Ke(t), (4) 

where e(t) = yi(t) — y2{t) denotes an output synchro- 
nization error and K is a scalar controller gain. The 
problem of finding static output feedback even for lin- 
ear systems is one of the classical problems of control 
theory. Although substantial effort has been devoted to 
its solution and various necessary and sufficient condi- 
tions for stabilizability by static output feedback have 
been obtained, most existing conditions are not testable 
practically [2J, [25|. In this paper we analyze a natu- 
ral and relatively broad class of systems for which con- 
structive conditions for output feedback stabilization are 
known is the class of passifiable (or feedback passive) sys- 
tems (for linear systems this introduced and studied in 
[22T [2a|). Since we are dealing with a nonlinear problem 
further complicated by information constraints, we re- 
strict our attention to sufficient conditions for solvability 
of the problem and evaluate upper bounds for synchro- 
nization error. 



III. CODING PROCEDURES 

In the paper [l9| the properties of observer-based syn- 
chronization for Lurie systems over a limited data rate 
communication channel with a one-step memory time- 
varying coder are studied. It is shown that an upper 
bound on the limit synchronization error is proportional 
to a certain upper bound on the transmission error. Un- 
der the assumption that a sampling time may be prop- 
erly chosen, optimality of binary coding in the sense of 
demanded transmission rate is established, and the rela- 
tionship between synchronization accuracy and an opti- 
mal sampling time is found. On the basis of these results, 
the present paper deals with a binary coding procedure. 

Consider the memoryless (static) binary quantizer to 
be a discretized map q : M. — > M as 

q(y,M) = Msign{y), (5) 

where sign(-) is the signum function: sign(y) = 1, if 
y > 0, sign(y) = — 1, if y < 0. Parameter M may be re- 
ferred to as the quantizer range. Notice that for a binary 
coder each codeword symbol contains one bit of informa- 
tion. The discretized output of the considered quantizer 
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is given as y = q(y,M). We assume that the coder and 
decoder make decisions based on the same information. 
The output signal of the quantizer is represented as a 
one-bit information symbol from the coding alphabet S 
and transmitted over the communication channel to the 
decoder. 

In time-varying quantizers 0, 0; [H, H3, HH the range 
M is updated with time and different values of M are 
used at each step, M = M[k]. Using such a "zooming" 
strategy it is possible to increase coder accuracy in the 
steady-state mode and at the same time, to prevent coder 
saturation at the beginning of the process [271 ] . 

In the present paper we use the following time-based 
zooming strategy for a quantizer range 



M[k} = M oP k , fc = 0,l,. 



(6) 



where < p < 1 is the decay parameter. The initial 
value Mo should be large enough to capture the region of 
possible initial values of y . Equations (J5J, ^ describe 
the coder algorithm. A similar algorithm is realized by 
the decoder. Namely, the sequence M[k] is reproduced 
at the receiver node utilizing © such that the values of 
y[k) are restored with the given M[k] using the received 
codeword s[k] S S. 



IV. EVALUATION OF SYNCHRONIZATION 
ERROR 

Let us evaluate the limit synchronization error, tak- 
ing into account transmission of the error signal over the 
communication channel and coding procedure. Since the 
control signal is piecewise constant over sampling inter- 
vals [tk,tk+i], the control law ([J| becomes 



u{t) = Ke{t), 



(J) 



where e(t) = e[k] as kT s < t < (k + 1)T S , e[k] is the 
result of transmission of the synchronization error signal 

e{t) = yi(t) — 2/2(0 over the channel, k = 0, 1, 

According to the quantization algorithm (J5]) , the quan- 
tized error signal e[k] becomes 



e[k]=M[k] sign(e(i fe )), 



(8) 



where e(t) = yi(t) — 2/2 (i), ife = kT s , and the range M[k] 
is defined by ©. 

Taking into account the stepwise shape of the control 
function in ([7]), rewrite the controller model in the fol- 
lowing form: 



u{t) = Ke(t)-KS(t), 



(9) 



where S(t) = S q (t) + 5 s (t) is a total error, S q (t) = e(tk) — 
e[k) = Ce(tfe)— e[k] is a quantization error, S s (t) = e(t) — 
e{tk) = Ceit) — Ce(tk) is a sampling error. 

It is seen from the quantization procedure ((5J) that if 
the value e(tk) satisfies the inequality |e(ifc)| < 2A/[fc], 
then the quantization error does not exceed M[k\: 



8 q (t)\ < M[k]. For the sampling error the following 
bounds hold: 



\5 s (t)\< / |£(r)|dr< max \i(r)\-T s . 

J t k <T<t 
tfc 

Therefore, the total error S(t) = 5 q (t) + S s (t) in the case 
when e(ifc) < 2M[k] is overbounded as: 

\5(t)\ <M[k}+ sup \e(s)\-(t-t k ). (10) 

t k <s<t 

A difficulty for evaluation of the synchronization error 
is in the dependence of right-hand side of (fTU|) on the 
trajectory of the system which is not known a priori. 
Therefore in order to evaluate the error dynamics we need 
to analyze dynamics of 6(t) and e(t) simultaneously. 

In order to analyze the synchronization error we make 
two assumptions: 

Al. Nonlinearity <p(y) is Lipschitz continuous: 

\f(Vi) ~ V>(V2)\ < L v \yi - 2/2 j (11) 

for all y%, 2/2 and some L v > 0. 

A2. The linear part of (pj is strictly passifiable: ac- 
cording to the Passification Theorem [13, [HI (see 
Appendix [A|) . this means that the numerator /3(A) 
of the transfer function W(X) = C(AI - A)^ 1 B 
= /3(A)/a(A) is Hurwitz (stable) polynomial of de- 
gree n — 1 with positive coefficients (the so-called 
hyper- minimum- phase (HMP) property). 

It follows from condition A2 and the Passification The- 
orem (see Appendix that the stability degree 770 of the 
polynomial /3(A) (a minimum distance from its roots to 
the imaginary axis) is positive and for any 7/: < r\ < t/q 
there exist a positive definite matrix P = P T > and a 
number K such that the following matrix relations hold: 

PA K + A T K P < -2r)P, PB = C t , A K = A- BKC. 

(12) 

Any sufficiently large real number can be chosen as the 
value of K. 

The main result of this Section is formulated as follows. 

Theorem 1. Let Al, A2 hold, the coder parameters T s , 
p and the auxiliary parameters 77, 77' be chosen in order 
to meet the inequalities 



T s b < 1, q<p<l, < 77' < 77 < 770, 



(13) 



where 

q = max I exp(— rj'T s ) - 

T s a (K + L v ) 
2(r ? -7 7 ')(l-TA)V% 



T s a Q (K + L v )(l - exp(-r)'T s )) 
2r/(l - T s b )VKZ 



4 



tjq is a stability degree of the polynomial /5(A). Let the 
number K is chosen from (|12[) and A m ; n is the minimum 
eigenvalue of the matrix P from (|T2|) . 

Let the coder range M[k] be chosen as follows 



M[k] = M oP k , 



(14) 



where 



Mn = W, 



0" 



( p- g)(l - T s b Q )y/X min 

r(i^ + L v ) 
1 - exp(- V 'T s ) 1 

2r/ ' 2(t) - r/) 



Then for all e(0) such that e(0) T Pe(0) < W the cur- 
rent values of the synchronization errors decrease expo- 
nentially: 



|#]| < ||e[*]|| < 2p k M . 



(15) 



Proof. The key point of the proof is comparison of the 
hybrid system in question with an auxiliary continuous- 
time system (the continuous model) possessing useful sta- 
bility and passivity properties (29l , |30| . 

Rewrite the error equation in the following form: 



e = Ae + B((e, t) - Bu, e = Ce, 



(16) 



where ((e,t) — cp(y(t)) —ip(y(t) — ej satisfies the inequal- 
ity |C(e,t)| <L v \e\. 

Substituting © into (jTHJ) we obtain 



e = A K e + BC(e, t) + BK8(t), 



(17) 



where Ax = A — BKC. Employing the HMP condition 
and the Passification Theorem, see Appendix [5] pick up 
the (n x n)-matrix P = P T > and the positive number 
K such that PA K + Aj { P < -2r}P, PB = C* T , and 
choose the Lyapunov function candidate V(e) = ^e T Pe. 
Introducing a new nonlinearity £(e,t) = C( £ j0 + L v e, 
satisfying the sector condition £e > and making the 
change K — > K + L v , transform equation (Tl7|) to the 
form 

e = A K e + B£-B(K + L v )S(t). (18) 

The time derivative of V(e) is evaluated as follows: 

V = e T (PA K + A T K P)e - e T PB(K + L v )S(t), 

or, after simple algebra 

V<-2nV + \e\{K + L v )\5{t)\, (19) 

where the equivalent disturbance S(t) satisfies the in- 
equality 



\6(t)\<M[k] 



sup 



\Ce(r)\\-T s . 



Overbounding ||e(r)|| from (| 18|) and taking into ac- 
count that ||e|| < y/V(e)/X m i n , where A m i n is the mini- 
mum eigenvalue of P yields 



\5{t)\ < M[k] + T s a JV[k + 1] + b S[k + 1 



(20) 



where V[k + 1] = sup V(e(t)), S[k + 1] = 
t k <t<t k+1 

sup \6(t)\,a = (\\CA K \\+\CB\-\\C\\)/V>^, 

b = \CB\ ■ (K + L v ). 

Let T s b < 1, i.e. T s < (\CB\ ■ (K + L V ))~ X . Taking 
the supremum of the left part of (|20[) over [tfc,tjfc+i] we 
arrive at an inequality 



M[k 



T s a 



Substituting (j2"Tj) into p^)) we obtain the functional- 
differential inequality for tk < t < tk+i 



V<-2r]V+Jv(e{t))(a sup y/V(t') + b[k]) , fc = 0, 1, . 
v t k <t'<t 

(22) 



where 



T s oq(K + L ¥ 

(i - &or fl )v% 



6[fc] 



M[fc](ff + L y ) 

(1 — &0^s)VAmin 



To estimate solutions of (|22|) employ Lemma 1 (see Ap- 
pendix [B|) with V(t) = V(e(t)Y Choosing the coder pa- 
rameters satisfying conditions of the Theorem and check- 
ing the conditions of Lemma 1 complete the proof. 

Based on Theorem 1 the following design method is 
proposed. The transmission rate T s should be chosen 
based on condition Al. Parameter p in (JHJ) should be 
found independently on T s based on the accessible sta- 
bility degree rj of the continuous model (|18|) (the max- 
imal value of r\ corresponds to the stability degree r]o 
of the numerator /3(A) of the transfer function W(X) = 
C(AI — A)~ 1 B. Finally, the value of r\ should be found 
based on the solution of LMI (fT2"|) for the chosen control 
gain K. 

Remark 1. Conditions q < 1, T s bo < 1 are always 
true for sufficiently small T s > (i.e. for sufficiently 
large capacity of the communication channel). Indeed, 
the following relations hold up to the second terms in T s : 



q w max < 1 — ?/T. 



T^a (K + L v ) T s a (K + L v ) 



2\^Xm 

The threshold for T, will be 



2(f) - V')VK 



T s < 



1 



K 



hn^r/v/A 



mm "^o ' 



2( V - V ')^/K~-a -\ \CB\- 1 }. 



(23) 
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It follows from ((23)) that grows of the controller gain 
K leads to increase of the transmission rate required for 
decay of the synchronization error. 

The dependence (|15|) will be used for numerical anal- 
ysis in Section fVl 

Remark 2. In stochastic framework the estimates of 
the mean square value of the synchronization error can 
be obtained. There is a significant body of work in which 
the quantization error signal S(t) is modeled as an extra 
additive white noise. This assumption, typical for digital 
filtering theory, is reasonable if the quantizer resolution 
is high [3l|, but it needs modification for the case of a 
low number of quantization levels • 

Remark 3. For practice, it is reasonable to choose the 
coder range M[k] separated from zero. The following 
zooming strategy for a quantizer range may be recom- 
mended instead of ([6]): 

M[k] = (M - A/ oo ) j0 fe + A/ oo , fc = 0,l,..., (24) 

where < M M < M stands for the limit value of M[k). 



V. EXAMPLE. SYNCHRONIZATION OF THE 
CHAOTIC CHUA SYSTEMS 

Let us apply the above results to synchronization of 
two chaotic Chua systems coupled via a channel with 
limited capacity. 

Master system. Let the master system (JTJ) be repre- 
sented by the following Chua system: 

Xx =p(-X! + if{yi) +x 2 ), t > 0, 
±2 = X\ - x 2 + x 3 (25) 
^3 = -qx 2 , 
yi (i) = xi{t), 

where y\(t) is the master system output, p, q are known 
parameters, x — [xi,X2,X3] T G R 3 is the state vector; 
tp(yi) is a piecewise-linear function, having the form: 



¥>(y) = m y + mifly + l\-\y- 1|), 



(26) 



where mo, mi are given parameters. 

Evidently, Chua system (|2"5")l may be represented in 
Lurie form ([T]) with the matrices: 



.4 



-p 


p 


o" 




P 


1 


-1 


l 


, B = 








-q 











C= [1,0,0]. (27) 



It is easy to check that the linear part of Chua system sat- 
isfies the HMP condition. Indeed, for the triple (A, B, C) 
from (EU) the transfer function W(X) = C(AI- A)- X B = 
/3(A)/a(A) is as follows: 



W(X) 



p{\ 2 + X + q) 



The numerator /3(A) = p(X 2 + X + q) is a Hurwitz poly- 
nomial of degree 2, i.e. the HMP condition holds for all 
p > 0, q > 0. 

Slave system. Correspondingly, the slave system equa- 
tions ([2]) for the considered case becomes 




t > 0, 



(28) 



where y2(t) is the slave system output, z = [z-\_, z%Y G 
K 3 is the state vector, ^(2/2) is defined by (|2"6"j) . 

Controller has a form ([7]), where the control gain K is 
a design parameter. 

Coding procedure has a form ([5]) , © . The input signal 
of the coder is e(t). The error signal e(t) of the controller 
([7|) is found by holding the value of e[k] over the sampling 

interval [kT S) (k + 1)T S ), k = 0, 1, The initial value 

Mq of the coder range and the decay factor p in ^ are 
design parameters. 



ti=0.3, R=25 bit/s 




A 3 + (1 +p)X 2 + qX+pq' 



FIG. 1: Time histories: X2(t) (dash line), Z2(t) (solid line) 
and synchronization error e^ii) = X2(t) — 22 (t) for n — 0.3, 
R = 25 bit/s. 

The following parameter values were chosen for simu- 
lation: 

- Chua system parameters: p = 10, q = 15.6, mo = 
0.33, mi = 0.945; 

- the controller gain K = 10. Feasibility of relations 
(|A2|) for this value of K and the given matrices A, 
B, C is checked by means of Yalmip package {32}; 

- the sampling time T s was taken from the interval 
T s G [0.02,0.1] s for different simulation runs (a 
corresponding interval for the transmission rate R 
is Re [10,50] bit/s); 
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- the coder parameters: Mo = 5, p = exp(—r]T s ), 
7] = 0.3; 

- the initial conditions for the master and slave sys- 
tems were: x = [3, — 1, 0.3] T , z — 0; 

- the simulation final time t^ n = 1000 s. 
The normalized state synchronization error 

max ||e(f)|| 
max \\x(t)\\ 

0<t<t tin 

where S y (t) = yi(t) — yi(t), e(t) = x{t) — z(t) was calcu- 
lated. 

Q 



10" - ; 

1CT 4 - | 

10" - i 



10 -"i ' — ^ < < < ' 

20 25 30 35 40 45 50 

R, bit/s 

FIG. 2: Normalized synchronization error Q v.s. transmission 
rate R. 

Simulation results are plotted in Figs. I1I2I 

Synchronization performance may be evaluated based 
on time histories of the state variables X2(t), zi(£) and 
the synchronization error e2(i). Typical trajectories are 
depicted in Fig. [TJ As seen from the plots, the synchro- 
nization transient time is about 15 seconds, which agrees 
with the chosen value of the coder parameter rj. 

The logarithmic graph of the normalized synchroniza- 
tion error Q as a function of the transmission rate R is 
shown in Fig. [21 It is seen from this plot that if the trans- 
mission rate exceeds the minimal bound R m m ~ 23 bit/s, 
the proposed controlled synchronization strategy ensures 
asymptotical vanishing the synchronization error. If the 
transmission rate is less that the bound R m i n , the syn- 
chronization is not always possible. 

Remark 4- An idealized problem has been considered 
in this paper to highlight the effect of the data-rate lim- 
itations in the closed-loop synchronization of nonlinear 
systems. In real- world problems external disturbances, 
measuring errors and channel imperfections should be 



taken into account. Evidently in the presence of irreg- 
ular nonvanishing, asymptotic convergence of the master 
and slave systems trajectories cannot be achieved. 

Remark 5. Similar results are obtained if the control 
signal is also subjected to information constraints. 

Remark 6. Optimality of the binary coder for synchro- 
nization under information constraints was established 
in 19] for the case when the master system output yi(t) 
rather than the output synchronization error e(t) is trans- 
mitted over the channel. The problem of coder optimiza- 
tion for the considered case is under investigation. 

Remark 7. The results obtained for synchronization 
problem can be extended to tracking problems in a 
straightforward manner, if the reference signal is de- 
scribed by an external (exogenious) state space model. 



VI. CONCLUSION 

Limit possibilities of controlled synchronization sys- 
tems under information constraints imposed by limited 
information capacity of the coupling channel are evalu- 
ated. It is shown that the framework proposed in 
is suitable not only for observer-based synchronization 
but also for controlled master-slave synchronization via 
a communication channel with limited information ca- 
pacity. 

We propose a simple coder-decoder scheme and pro- 
vide theoretical analysis for multi-dimensional master- 
slave systems represented in Lurie form. An output feed- 
back control law is proposed based on the Passification 
Theorem [22|, [23|. It is shown that the synchronization er- 
ror exponentially tends to zero for sufficiently high trans- 
mission rate (channel capacity) . The key point of the syn- 
chronization analysis is comparison of the hybrid system 
in question with an auxiliary continuous-time system (the 
continuous model) possessing useful stability and passiv- 
ity properties. Such an approach was systematically de- 
veloped in the 1970s under the name of the Method of 
Continuous Models (29l . [30j . 

The results are applied to controlled synchronization of 
two chaotic Chua systems via a communication channel 
with limited capacity. Simulation results illustrate and 
confirm the theoretical analysis. 

Unlike many known papers on control of nonlinear 
systems over a limited-band communication channel, we 
propos and justify a simple coder/decoder scheme, which 
does not require transmission of the full system state vec- 
tor over the channel. A constructive design method for 
controller and coder/decoder pair is proposed and esti- 
mates of the convergence rate are given. The results ob- 
tained for synchronization problem can be extended to 
tracking problems in a straightforward manner, if the 
reference signal is described by an external (exogenious) 
state space model. 

Future research is aimed at examination of more com- 
plex system configurations, where channel imperfections 
(drops, errors, delays) will be taken into account. 
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APPENDIX A: PASSIFICATION THEOREM 



Then the following statements hold: 



Consider a linear system 



e = Ae + B£(t), e = Ce 



(Al) 



with transfer function W(X) = C(XI - A)~ 1 B= 
6(A)/a(A), where b(X), a(X) are polynomials, degree of 
a(A) is n, degree of 6(A) is not greater than n — 1. The 
system is called hyper-minimum phase (HMP), if 6(A) 
is Hurwitz polynomial of degree n — 1 with positive co- 
efficients. To find existence conditions for a quadratic 
Lyapunov function we need the following result. 

Passification Theorem |H, [jil /. There exist positive- 
definite matrix P = P T > and a number K such that 

PA K + A T K P<0, PB = C T , A K = A - BK C (A2) 

if and only if W(X) is HMP. 

Consider a linear system with feedback 

e = A K e + B£(t), e = Ce, A K = A- BKC. (A3) 

Remark Al. It follows from the Passification Theorem 
that there exist a quadratic form V(e) = e T Pe and a 
number K such that time derivative V(e) of V(e) along 
trajectories of (|A3|) satisfies relation 



V(e) < for & > 0, x ^ 



(A4) 



if and only if M^A) is HMP. Indeed, assume that K is 
fixed. Relation (IA4[) is equivalent to existence of the 
matrix P = P T > such that e T P{A K e + BS,) +£Ce < 
for x ^ 0. Since £ is arbitrary, the latter in turn, is 
equivalent to matrix relations PAk + A\P < 0, PB = 
C T and, by Passification Theorem, to HMP condition. 

Remark A2. It also follows from Passification Theorem 
that if HMP condition holds then K satisfying (|A2[) can 
be chosen sufficiently large. Besides, zero matrix in the 
right hand side of the inequality in (|A2|) can be replaced 
by matrix — jj,P for sufficiently small /j > 0. 



APPENDIX B: LEMMA 1 



Lemma 1. Consider the functional-differential inequal- 



ity 



V(t) < -2t]V(t) + VW)(a sup ^/VW) + b[k]) , 

V t k <f<t J 

t k <t<t k+1 , t k = kT s , k — 0,1, ... . (Bl) 
Let < r/ and q < 1, where 

a(l - exp(-?/T s )) a } 



q = max < exp(— f/T s ) + 



2r/ 



2(r) - if) j ' 
(B2) 



1. W[k + l]<qW[k]+rb[k], 

1 - exp(-r/T s ) 1 



where r = max , 

2rj' 

fc-i 



2(77-7/)]' 



(B3) 



(B4) 



i=0 



3. If q < p < 1 and b[k] = 



p-q 



W[k] 



then W[k] decays exponentially: W[k] — p k Wo- 
Proof of Lemma 1. 



Make the change W(t) = \/V(t) and introduce the 

notation W[k] = W{t k ), W(t) = max W(t). Then 

tk-\<t<t 

inequality (|Bip becomes 

n h\h] 

(B5) 



W < -r/W+^W(t) + 



Let us prove an auxiliary statement: for any ?/: < 
i]' < 77 and for any k = 0, 1, . . . the function W{t) de- 
creases and the following inequality holds: 

W(t) < max I.W{T k ) exp ( - V '(t - t k )) 



+ 



aW{t ^ m (l-eM-V(t-t k ))), 



aW{t k ) + b[k] 
2(»7 ~ rf) 

To this end consider the inequality 



(B6) 



-VW(t) + \W{£) + M < ^^(f). (B7) 

It holds for t = t k , and, obviously holds for some in- 
terval t k <<<<'. Let if be the maximum of instants for 
which dHZj) holds. It follows from (JB5J) that W(t) < 
and VT(t) < W{t k ) for t k < t < t'. Hence, W(t) = W(t) 
for t k < t < t' , and (|T37|) becomes 



f^ fe ) + M< (7? _ VW) 



and it holds iff 



W(t) > 



aW{t k )+b[k] 
2(7? - »?') 



(B8) 



(B9) 



Therefore, as long as (|B9[) holds, the value of W(t) 
decreases exponentially and 

W(t)< W(t k )exp(-rf(t-t k )) 
aW(t k )+b[k] 



+ 



2rf 



1 - exp ( - 7/(t - tfc)) . (B10) 



On the other hand, if the opposite inequality W(t) < 
aW 2( v -nT ] holds for some t = t '> [t holds for a11 t:t' < 
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t < tk+i owing to monotone decrease of W(t). Therefore, 
(|B6I) is proven. 

The first statement of the lemma follows from (|B6j) if 
we set t = tk+i, collect terms with W(tk) and b[k] and 
take their maximum values. The second statement is 
derived by iteration of (|B3|) and the third one is obtained 
by direct substitution. 
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